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FOREWORD 

This  report  describes  an  analytical  investigation  of  the  applicability 
of  surface  tension  screens  for  use  in  ramjet  fuel  systems  using  the  heavy 
hydrocarbon  fuel  RJ-5.  The  analysis  was  performed  in  the  Air  Force  Aero 
Propulsion  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio  by  Jack  R. 
Fultz  (RJA).  The  work  was  performed  under  Project  3012,  "Ramjet 
Technology",  Task  301211,  "Ramjet  Design  and  Assessment",  Work  Unit 
Number  30121103,  "Multi-Purpose  Missile  Design  Studies"  during  the 
period  January  1973  through  July  1973. 

This  report  was  submitted  by  the  author  April  1974. 

This  technical  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 

This  report  presents  an  analysis  of  surface  tension  screens  for 
their  applicability  to  fuel  systems  in  ramjet-powered  vehicles.  Selection 
of  surface  tension  screens  for  fuel  expulsion  over  other  expulsion 
techniques  was  made  on  the  basis  of  their  high  temperature  capability, 
moderate  pressure  requirements,  geometrical  considerations,  simplicity 
(no  moving  parts),  and  other  factors.  Analysis  of  the  pressure  losses 
encountered  in  expelling  the  heavy  hydrocarbon  fuel  (RJ-5)  used  in  this 
study  revealed  that  the  primary  pressure  loss  in  the  system  was  encountered 
in  flowing  the  fuel  through  the  surface  tension  screen.  The  time  required 
to  expell  fuel  from  various  trap  tank  sizes  was  calculated  as  a function 
of  fuel  flow  rate.  Parameters  investigated  in  addition  to  trap  tank  size 
included  orifice  vent  area  and  acceleration  (g's)  force.  Result  of  this 
analysis  revealed  that  large  pressure  drops  were  encountered  in  flowing 
the  RJ-5  fuel  through  surface  tension  screens  having  sufficiently  small 
holes  to  provide  appreciable  surface  tension  force.  The  primary  reason 
for  this  large  pressure  drop  is  the  high  viscosity  of  the  RO-5  fuel. 

The  utility  of  surface  tension  screens  for  ramjet  fuel  systems  using 
RJ-5  fuel  is  limited  to  high  fuel  tempera^ires  at  low  (^3)  acceleration 
(g's)  level.  If  low  viscosity  blends  of  high  density  fuels  are  developed 
and  further  analysis  reveals  that  low  g operation  is  feasible  during  the 
critical  fuel  expulsion  cycle,  then  a more  detailed  analysis  of  surface 
tension  screens  for  ramjet-powered  vehicles  will  be  warranted. 
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SECTION  I 
INTRODUCTION 

1.  GENERAL  CONSIDERATIONS  FOR  SELECTING  SURFACE  TENSION  SCREENS 

Fuel  expulsion  from  the  fuel  tanks  of  ramjet-powered  missiles  has 
generally  been  accomplished  by  positive  expulsion  techniques,  primarily 
elastromeric  bladders  or  diaphragms  pressurized  by  stored  gas  bottles  or 
ram  air.  The  missiles  which  have  employed  this  technique  have  operated 
over  moderate  flight  envelopes,  primarily  below  Mach  3.  In  this  flight 
envelope  the  aerodynamic  heating  and  the  ram  air  recovery  temperatures 
were  low  enough  to  permit  use  of  elastromeric  bladders.  Overall  expulsion 
efficiencies  achieved  with  these  elastromeric  systems  were  95%  or  better. 
Although  elastromeric  expulsion  techniques  have  proven  useful  in  a number 
of  ramjet-powered  missiles  their  utility  becomes  less  apparent  for  future 
systems  operating  at  Mach  numbers  considerably  above  Mach  3.  The  thermal 
load  imposed  on  the  bladder  above  Mach  3 becomes  so  great  that  the 
structural  integrity  of  the  material  becomes  a problem.  An  obvious  method 
of  overcoming  the  thermal  problem  is  to  use  a layer  of  insulation  between 
the  elastromeric  material  and  the  metal  skin  of  the  missile.  This  imposes 
a weight  and  volume  penalty  for  a given  missile  configuration  resulting 
in  some  decrease  in  performance  of  the  missile. 

Other  expulsion  techniques,  shown  in  Figure  1,  such  as  bellows, 
metallic  bladders,  pistons,  and  surface  tension  screens  have  been  used 
for  orientation  and  expulsion  of  fuel  and  oxidizers  in  a number  of  missile 
and  space  vehicle  propulsion  systems.  All  of  these  expulsion  devices  are 
driven  by  a pressurizing  gas  such  as  air  or  nitrogen  (Reference  1).  A 
cursory  examination  of  these  alternate  expulsion  techniques  was  performed 
for  their  application  to  the  multipurpose  missile  (MPM).  Although  several 
of  these  expulsion  techniques  may  be  attractive,  various  factors  such  as 
pressure  requirements  and  tank  shape  resulted  !n  the  selection  of  the 
surface  tension  screen  concept  for  further  in-house  study. 
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Figure  Id  is  a schematic  of  a typical  surface  tension  screen  system. 

One  or  more  screens  having  hole  sizes  in  the  micronic  range  are  contained 
inside  a trap  tank  which  holds  fuel  in  contract  with  the  screens  during 
vehicle  maneuvers  which  would  otherwise  rapidly  uncover  the  surface 
tension  screens.  Fuel  flows  through  the  screen(s)  by  means  of  pressurizing 
the  fuel  tank  with  an  inert  gas  such  as  nitrogen.  The  screen(s)  permit 
only  liquid  to  flow  until  the  total  pressure  drop  across  the  screen 
orifices  exceeds  the  surface  tension  force  of  the  liquid  under  consideration. 
(See  Section  II  for  operation  principle).  If  this  surface  tension  force 
is  exceeded,  gas  ingestion  in  the  screen (s)  will  occur  ana  two-phase  flow 
to  the  engine  will  result.  The  trap  tank  containing  the  surface  tension 
screen(s)  is  perforated  to  allow  fuel  flow  into  the  screen(s)  and  to 
permit  the  pressurizing  gas  to  escape  when  the  vehicle  returns  to  normal 
"G"  operation. 

An  analytical  evaluation  of  the  surface  tension  screen  concept  for 
application  to  future  ramjet- powered  missiles  has  been  performed.  Primary 
advantages  of  the  surface  tension  screen  over  the  bladder  expulsion  system 
are: 


a.  less  inert  weight  and  volume 

b.  no  moving  parts 

c.  no  thermal  limitation  (same  as  total  missile)  for  the  device  itself* 

d.  unlimited  expulsion  cycle  capability  (Reference  2). 

2.  DESIGN  CONSIDERATIONS  FOR  SELECTION  OF  SURFACE  TENSION  SCREENS 

The  initial  (and  still  primary)  use  of  surface  tension  devices  was 
for  space  vehicles  under  conditions  of  low  or  near  zero  gravity.  The 
three  primary  forces  which  influence  the  dynamic  behavior  of  liquid/gas 
systems  are  body  forces,  capillary  forces,  and  viscous  forces.  The 
surface  tension  screen  concept  is  effective  only  in  systems  where  the 
capillary  forces  predominate.  The  relative  importance  of  the  body  and 

*Temperature  will  affect  fluid  properties  which  will  affect  efficiency 
of  the  surface  tension  device. 
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capillary  forces  for  a specific  case  (i.e.,  evaluation  of  surface  tension 
devices)  can  be  estimated  by  three  interrelated  dimensionless  parameters 
known  as  the  Weber  number,  the  Froude  number,  and  the  Bond  number.  The 
Weber  number  is  defined  as  the  ratio  of  inertial  force  to  capillary  force: 


W 


e 


0) 


where,  W = 
P = 


Weber  number,  dimensionless 

3 

liquid  density,  lb-m/ft 


V = liquid  velocity  through  the  surface  tension  screen,  ft/sec 
d = hole  diameter  of  the  surface  tension  screen,  ft 


o = surface  tension  of  the  liquid,  Ib-f/ft 

g = gravitational  constant  = 32.174  Ik^ILLi 
c lb-f  sec' 


For  Weber  numbers  much  greater  than  1,  capillary  forces  are  insignificant 
so  that  liquid  motion  in  the  system  is  determined  by  the  inertia  of  the 
system.  For  Weber  numbers  much  less  than  1,  capillary  forces  dominate 
and  define  liquid  movement  in  the  system. 


The  Froude  number  provides  an  estimate  of  the  inertia  force  to  the 


gravitational  force  and  is  defined  as  follows: 


(2) 


where  F^  = 
V = 

9 B 
d = 


Froude  Number,  dimensionless 

liquid  velocity  through  the  surface  tension  screen,  ft/sec 

2 

local  acceleration,  ft/sec 
hole  diameter,  ft 


For  Froude  numbers  very  much  greater  than  1.0  the  gravitational  forces 
are  sufficiently  weak  that  they  have  essentially  no  influence  on  the 
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fluid  motion  whereas  for  very  small  values  of  the  Froude  number  the 
gravitational  forces  must  be  considered. 

The  Bond  number  is  the  ratio  of  the  gravitational  forces  to  the 
capillary  forces  and  is  obtained  by  dividing  the  Weber  number  by  the  Froude 
number. 


B 


o 


= _ pgd£ 


gd 


(3) 


For  very  small  values  of  the  Bond  number,  the  capillary  forces  predominate 

in  determing  liquid  motion  whereas  the  gravitational  force  predominates 

for  large  Bond  numbers.  Figure  2 depicts  the  three  hydrodynamic  regimes 

defined  by  the  Froude,  Weber,  and  Bond  numbers.  A number  of  calculations 

were  performed  to  evaluate  the-  hydrodynamic  regime  anticipated  for  the 

surface  tension  screens  considered  by  this  analysis.  These  calculations 

were  performed  for  a single  screen  element.  The  screen  considered  v/as 

one  currently  available  from  the  Western  Filter  Company  and  was  evaluated 

by  Vought  Missile  and  Space  Company  (VMSC)  for  use  in  the  Harpoon  missile. 

The  total  screen  area  for  this  particular  element  is  160  in^  (1.1;  I ft^) 

2 

with  a porosity  of  30.7%.  The  fuel  flow  area  therefore  is  0.3411  ft  . 

The  screen  tested  by  VMSC  was  a 10  micron  unit.  For  this  analysis,  44 

and  100  micron  screen  sizes  were  assumed  in  addition  to  the  10  micron 

2 

baseline  screen.  Flow  area  was  kept  constant  at  0.3411  ft  for  the 
larger  hole  sizes.  Bond  numbers  and  Weber  numbers  were  calculated  for 
the  three  hole  sizes  for  fuel  flow  rates  of  1.0  and  10  pounds  per  second 
and  for  acceleration  forces  of  1.0  and  10  g's.  Results  are  presented  in 
Table  I. 


-5  -2 

The  Bond  numbers  obtained  ranged  from  2.64  x 10  up  to  2.45  x 10  . 

-5  -2 

Weber  numbers  obtained  ranged  from  2.77  x 10  up  to  4.85  x 10  . 
Referring  to  Figure  2,  it  is  seen  that  these  low  values  of  the  Bond 
and  Weber  numbers  place  the  flow  through  these  surface  tension  screens 
in  the  capillary-dominated  regime.  The  primary  reason  for  these  low 
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values  of  the  Bond  and  Weber  numbers  in  the  surface  tension  screens  is 
the  extremely  small  hole  diameters  considered.  This  is  seen  by  examining 
Equation  1 for  the  Weber  number  where  the  hole  diameter  is  to  the  first 
pcwer  in  the  numerator  and  by  examining  Equation  3 for  the  Bond  number 
where  the  hole  diameter  is  to  the  second  power  in  the  numerator.  Also 
contributing  to  the  low  Weber  numbers  is  the  fact  that  the  velocity 
through  the  screens  is  less  than  1 ft/sac.  Since  the  velocity  is  to  the 
second  power  in  the  numerator,  the  velocity  term  becomes  extremely 
important  in  producing  small  Weber  numbers  as  the  velocity  through  the 
screen  is  decreased  below  1 ft/sec. 

The  other  factors  used  to  define  these  dimensionless  parameters 
(Bond  and  Weber  numbers)  are  less  important  than  the  screen  dimensions. 
This  can  be  seen  by  examining  the  results  presented  in  Table  I.  Reducing 
the  fuel  temperature  from  77°F  to  0°F  produces  only  minor  changes  in 
the  Bond  and  Weber  numbers  as  does  increasing  acceleration  level  from 
1 g to  10  g's.  Thus,  one  could  conclude  that  these  screens  should 
perform  well  with  RJ-5  fuel  since  the  Bond  and  Weber  numbers  indicate 
that  the  flow  is  in  the  capillary-dominated  regime.  However,  one 
extremely  important  fuel  property  not  included  in  the  Bond  and  Weber 
numbers  is  viscosity.  As  will  be  shown  later,  this  property  is  critical 
to  the  performance  of  surface  tension  screens. 
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SECTION  II 

PRIMARY  DESIGN  PARAMETERS 
1 . PRINCIPAL  OF  OPERATION 

The  operation  of  surface  tension  screens  utilizes  capillary  force  as 
the  governing  parameter.  The  technique  generally  used  to  expel 1 fuel 
from  fuel  tanks  containing  surface  tension  screens  is  to  pressurize  the 
tank  with  an  inert  gas  to  a sufficiently  high  pressure  to  maintain  the 
fuel  flow  rate  required  by  the  engine.  Near  the  end  of  a flight  when  the 
fuel  supply  is  low  or  during  certain  maneuvers  of  the  missile  which 
produce  negative  G operation,  the  surface  tension  screen(s)  may  be 
exposed  to  the  pressurizing  gas.  Since  two  phase  flow  to  the  engine 
cannot  be  tolerated,  the  fuel  expulsion  system  must  be  designed  to  prevent 
gas  ingestion.  The  operation  of  a surface  tension  screen  is  shown 
schematically  in  Figure  3.  The  situation  depicted  in  Figure  3 is  for  a 
positive  g low  fuel  level  condition  which  results  in  partial  exposure  of 
the  surface  tension  screen.  The  principle  of  operation,  however,  is  not 
altered  for  negative  g orientation  of  the  fuel  tank. 

The  basic  condition  which  must  be  met  to  permit  gas-free  liquid  to 
flow  is  that  the  total  pressure  differential  between  the  gas  and  liquid 
at  all  points  across  the  screen  must  be  less  than  the  surface  tension 
pressure  force.  This  inequality  can  be  expressed  as  follows: 

APt  <Pcr  (4) 

where  AP^  = total  pressure  differential  across  surface  tension  screen 
orifice,  lb-f/ft^ 

APa  = surface  tension  pressure  drop  across  screen  orifice, 

1 b-f 
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Figure  3.  Operation  Principle  of  Surface  Tension  Screens 
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From  Figure  3,  it  is  seen  that  the  total  pressure  differential  across  a 
surface  tension  screen  orifice,  APy,  is  given  by  the  difference  in  the 
gas  pressure  outside  the  screen  minus  the  total  pressure  inside  the 
screen.  Expressed  mathematically,  this  pressure  differential  is  as 
follows: 


4Pr  ■ PG  - PTi 

where  P^.  is  the  total  internal  pressure. 


(5) 


If  the  gas  pressure,  Pg,  is  constant  then  the  pressure  differential  across 
the  screen  is  maximum  at  the  point  of  minimum  pressure  within  the  surface 
tension  screen.  From  Figure  3,  this  minimum  point  is  seen  to  occur  at 
the  outlet  plane,  P^,  since  the  maximum  flow  and  head  losses  have  occurred 
at  that  point.  Therefore,  Equation  (5)  becomes: 


- P 


3 


Equation  6 holds  as  long  as  two  phase  flow  does  not  occur. 


(6) 


The  pressure  at  P^  is  equal  to  the  gas  pressure,  Pg,  minus  the 
pressure  losses  which  occur  along  the  liquid  flow  path  (See  Figure  3). 
These  pressure  losses  are  determined  as  follows: 


P 


(7) 


where  P,  = total  pressure  at  any  given  point  in  trap  tank  outside  surface 

2 

tension  screen,  Ib-f/ft 

o 

Pg  - pressurant  gas  pressure,  lb-f/ft 

p 

h,  = pressure  of  liquid  head  at  any  given  point  in  trap  tank 

? 

outside  of  surface  tension  screen,  lb-f/ft 
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The  pressure  at  the  same  liquid  level  inside  the  surface  tension  screen 
is: 


P2  " P1  " APs 


P2  = PG  + % ' APs 


(8) 

(9) 


where  AP$  is  liquid  pressure  drop  across  the  surface  tension  screen 
orifice,  lb-f/ft2. 


The  pressure  at  the  screen  outlet  plane,  P^,  is 


P 


3 


or 


P3 


PG  + Ph, 


which  simplifies  to: 

P3  = PG  " Ph  ‘ *Ps  " APF 


(10) 


(ID 


where  Ph  = pressure  of  liquid  head  inside  screen  above  traj  tank  level, 
lb-f 
ft2 

2 

APp  = flow  friction  loss,  lb-f/ft 

Therefore,  from  Equation  6,  the  largest  pressure  differential  across 
the  portion  of  the  surface  tension  screen  exposed  to  gas  is: 

APt  = Pr  - (Pfi  - Ph  - AP  - APF)  (12) 

TMAX  G 6 h s F 

which  reduces  to 

APt  = Ph  + AP  + APp  /13) 

TMAX  h s F UJj 
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As  the  fuel  level  in  the  trap  tank  decreases,  each  of  the  terms  in 
Equation  13  increases,  assuming  positive  g operation  and  a constant  flow 
rate  out  of  the  surface  tension  screen.  That  this  is  true  can  be  seen  by 
the  following  reasoning: 

a.  Liquid  height  used  to  obtain  head  pressure  (P^)  increases  as  trap 
tank  level  decreases  resulting  in  higher  head  pressure. 

b.  For  a constant  flow  rate,  the  screen  flow  pressure  drop  increases 
with  decreasing  trap  tank  liquid  level  since  the  flow  area  decreases  with 
decreasing  liquid  level  of  the  trap  tank. 

c.  Flow  friction  pressure  drop  increases  with  decreasing  trap  tank 
liquid  height  since  the  percentage  of  fuel  which  must  travel  the  longer 
distance  to  the  outlet  is  increased. 


The  surface  tension  pressure,  AP^,  can  be  calculated  for  a given 
surface  tension  screen  by  the  following  equation  (Reference  3): 


*Po 


(14) 


2 

where  o = surface  tension  of  liquid,  Ib-f/ft 

d = hole  diameter  of  surface  tension  screen,  ft 


This  pressure  drop  can  be  expressed  in  terms  of  liquid  head  by  utilizing 
the  following  relationship: 


. h„  * 


(15) 


1 b-m 

where  p = density  of  liquid,  — £- 

ft6 

ha  - height  of  liquid,  ft 

2 

g = local  acceleration,  ft/sec 

g = gravitational  constant,  32.174  — , 

c Ib-f  sec 
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Rearranging  Equation  15  gives. 


h 


a 


Substituting  Equation  14  into  16  produces: 


(16) 


h _ 4a 

«'dl,  O?) 

Sc 

Equation  17  gives  the  maximum  pressure  in  terms  of  liquid  head  which  a 
surface  tension  screen  can  support  as  a function  of  liquid  properties 
(surface  tension  and  density),  hole  size,  and  acceleration  (g)  level, 
before  gas  flow  will  occur  across  the  screen.  Figure  4 shows  the  retention 
capability  of  various  screen  sizes  as  a function  of  temperature  and 
imposed  acceleration  levels  for  RJ-5  fuel.  From  Equation  17  it  is  seen 
that  retention  capability  is  inversely  proportional  to  hole  diameter  and 
to  acceleration  (g)  level.  For  example,  from  Figure  4,  the  liquid  head 
for  RJ-5  with  the  10  micron  size  screen  under  a 1 g load  is  rapidly 
decreased  to  a low  value  by  increasing  the  hole  diameter  or  by  increasing 
the  g loading.  Temperature  of  the  fuel  has  a less  drastic  effect  since 
the  decrease  in  retention  capability  from  -65°F  to  +500°F  is  only  about 
m of  the  -65°F  value. 


Equation  4 is  the  basic  inequality  which  must  exist  for  surface 
tension  screen  operation.  This  inequality  can  be  expressed  in  terms  of 
liquid  head  by  converting  the  pressure  drop  terms  to  liquid  head  terms. 
The  resulting  expression  is: 


hh  * hs  +hF  <ho 


(18) 


where  hu  = liquid  head  inside  screen  above  trap  tank  level,  ft 

hs  = equivalent  liquid  head  due  to  flow  through  screen,  ft 

hp  = equivalent  liquid  head  due  to  flow  friction  losses  inside 
due  to  surface  tension  of  the  fluid,  ft 
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TEMPERATURE  °F 

Figure  4.  Retention  Capability  of  Surface  Tension  Screens 
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Thus,  it  is  seen  that  in  the  operation  of  a surface  tension  screen  the 
total  liquid  head  due  to  height  and  flow  losses  must  be  less  than  the 
maximum  liquid  head  which  a screen  can  support  based  on  the  surface 
tension  of  the  liquid. 


The  most  difficult  factor  to  evaluate  is  the  equivalent  head  due  to 
pressure  drop  across  the  screen,  h$.  A number  of  papers  (References  4, 

5,  and  6)  have  been  published  giving  correlations  for  obtaining  screen 
pressure  drops.  Assumptions  used  in  these  correlations  have  included  the 
following  models  for  flow  through  the  screen: 

(1)  Orifices  or  Nozzles  in  Parallel  (Reference  4) 

The  flow  through  a screen  can  be  considered  as  flew  through 
a number  of  orifices  or  nozzles  in  parallel.  Thus,  the  pressure  drop  or 
head  loss  across  a screen  can  be  computed  from  an  orifice-type  equation. 
The  resulting  equation  for  head  loss  is: 


h 


s 


(19) 


where  h$  = head  loss  across  screen,  feet  of  fluid  flowing 

n = number  of  screens  in  series,  dimensionless 

C = screen  discharge  coefficient,  dimensionless 

« = fractional  free  projected  area  of  screen,  dimensionless 

V = superficial  velocity  ahead  of  screen,  ft./sec* 

2 

gc  = gravitational  acceleration,  32.17  ft/sec 


Experimental  data  [Volokhov,  Vestnik,  Ing.  Techn.,  (Reference  4), 
149-152  (1930)]  indicate  that  for  a series  of  screens  the  overall  head 
loss  is  directly  proportional  to  the  number  of  screens  in  series,  as 
given  by  Equation  19,  and  is  not  affected  by  either  the  spacing  between 
successive  screens  or  by  their  orientation  with  respect  to  one  another. 


where  Q = volumetric  flow  rate, 


ft° 

sec 


and  A = flow  area,  ft c 
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Screen  discharge  coefficient  "C"  is  a function  of  screen  Reynolds 
number. 


DoVp 

N = - 
NRe 


(20) 


where  = aperture  width,  ft 

p = fluid  density,  lb-m/ft3 
p = fluid  viscosity,  lb-f/(ft) (sec) 

For  plain  rectangular-mesh  screens,  Lapple's  plot  of  C vs  NRe  is  given 
in  Figure  :5.  This  curve  represents  most  of  the  data  to  within  20  percent. 
Coefficients  greater  than  1 probably  indicate  that  the  effective  free  area 
is  larger  than  that  of  the  projected  area  and  that  there  is  partial 
recovery  of  head  due  to  the  downstream  rounding  of  the  wires. 

(2)  Packed  Bed  (Reference  5) 

This  model  was  developed  for  flow  of  Newtonian  fluids  through 
all  types  of  woven  metal  screens.  The  screen  is  treat. J as  a very  thin 
packed  bed  in  which  the  pressure  drop  across  the  screen  is  considered  to 
be  the  sum  of  both  viscous  and  inertial  forces.  The  model  correlates 
friction  factor  versus  Reynolds  Number.  Pressure  drop  is  included  in  the 
friction  factor  term  so  that  once  the  friction  factor  is  obtained  from  a 
plot  of  Reynolds  Number  vs  friction  factor,  the  pressure  drop  through  the 
screen  can  be  calculated.  The  model  was  developed  by  using  experimental 
data  of  pressure  drop  vs.  velocity  for  nitrogen  over  a velocity  range 
from  0.1  to  30  ft/sec.  This  correlation  was  shown  to  hold  for  Reynolds 
Numbers  up  to  1000  for  water  as  the  test  fluid  by  substituting  data 
obtained  by  an  earlier  investigator  into  the  correlation  equation.  Results 
obtained  were  within  the  experimental  error  of  the  nitrogen  data. 
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Figure  5.  Screen  Discharge  Coefficients  for  Surface  Tension  Screens 
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(3) .  Win  Tec  Pressure  Drop  Model 

The  model  currently  being  used  by  McDonnell  Douglas 
Astronautics  Company  for  their  surface  tension  screen  work  on  the  Harpoon 
missile  was  developed  by  the  Win  Tec  Corporation  under  a NASA  contract 
(Reference  7).  This  model  correlated  screen  pressure  drop  data  for  JP-4, 
water,  hydrazine,  ethylene  glycol/water  mixtures,  liquid  nitrogen,  and 
hydraulic  oil  much  better  than  the  original  model  (Reference  8)  used  by 
McDonnell  Douglas. 

(4) .  Other  Methods 

Other  methods  of  correlating  overall  frictional  losses 
across  screens  are  given  by  Cornell  [Trans.  Am.  Soc.  Mech.  Engrs.,  80, 
791-799  (1958)]  and  Grootenhuis  [Proc.  Inst.  Mech.  Engrs.  168A,  837-846 
(1954)]. 

2.  SCREEN  PRESSURE  DROP  CALCULATIONS 

Pressure  drop  through  surface  tension  screens  as  a function  of  mass 
flow  rate  were  calculated  for  RJ-5  at  three  fuel  temperatures  using  the 
Reference  4 technique  described  in  subparagraph  1 above*.  Two  screen 
sizes  (44  and  100  microns)  were  used  and  all  calculations  were  for  the 
one  g level.  Results  are  shown  in  Figure  6.  The  following  example  will 
demonstrate  the  impact  of  the  Figure  6 data: 

EXAMPLE  1:  Use  44  micron  screen  at  1 g 
Fuel  Temperature  = 77°F 
From  Figure  4,  h , = 13.3  inches  of  RJ-5 

Equation  18  then  becomes:  13.3  > h^  + h$  + hp  or  13.1  - h$  > h^  + 

hp.  h$  is  obtained  from  Figure  6 assuming  maximum  flow  rate  of  10  Ib/sec. 
h$  = 1.05  ft  of  RJ-5  = 12.6  inches  of  RO-5.  Then,  13.3  - 12.6  > h^  + 
hp.  Thus,  the  total  allowable  pressure  drop  for  the  head  and  the  channel 


*The  Reference  4 technique  is  used  since  the  accuracy  required  in  this 
analysis  did  not  warrant  use  of  the  more  complex  technique. 
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flow  pressure  drops  is  only  0.7  of  an  inch  of  RJ-5.  Thus,  even  if  the 
channel  flow  pressure  drop  is  ignored  (a  reasonable  assumption)  the  total 
liquid  head  which  the  surface  tension  screen  can  withstand  is  only  0.7  of 
an  inch  of  RJ-5. 

Using  the  same  conditions  as  Example  1 except  for  fuel  temperature 
= 0°F,  it  is  seen  that  the  maximum  allowable  pressure  drop  is  exceeded: 

From  Figure  4,  ho  = 14.5  inches  of  RJ-5 

From  Figure  6,  h$  = 11.5  ft  of  RJ-5 

= 138  inches  of  RJ-5 

Then,  from  Equation  18,  h^  + 138  + hp  < 14.5. 

Thus,  the  allowable  pressure  drop  of  14.5  inches  of  RJ-5  is  overwhelmed 
by  the  138  inches  of  pressure  drop  through  the  screen.  This  means,  of 
course,  that  if  the  system  were  pressurized  to  overcome  this  screen 
pressure  drop,  ingestior,  of  the  pressurant  gas  would  occur. 

Even  though  the  pressure  drop  through  the  screen  decreases  with  fuel 
flow  rate,  the  screen  pressure  drop  still  exceeds  the  allowable  limit 
for  the  44  micron  screen  with  0°F  fuel  for  the  entire  range  of  fuel  flow 
rates  of  interest. 

The  primary  reason  for  the  excessive  pressure  drop  through  the  screen 
is  the  high  viscosity  of  the  RJ-5  fuel.  This  can  be  seen  by  examining 
Equation  18  which  shows  that  the  head  loss  through  the  screen  is  inversely 
proportional  to  the  square  of  the  screen  discharge  coefficient.  As 
shown  by  Figure  5,  the  screen  discharge  coefficient  is  a function  of  the 
Reynolds  Number  of  the  fuel  in  the  screen.  Because  Reynolds  Number  is 
inversely  proportional  to  the  fuel  viscosity  (see  Equation  19)  large 
values  of  viscosity  produce  small  Reynolds  Numbers.  From  Figure  5,  it 
is  seen  that  small  values  of  the  Reynolds  Number  correspond  to  small 
values  of  the  screen  discharge  coefficient  which,  in  turn,  produce  large 
pressure  drops  when  inserted  into  Equation  18.  For  this  study,  Reynolds 
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Numbers  through  the  screens  of  less  than  1.0  were  generally  obtained. 
These  low  values  of  the  Reynolds  Number  correspond  to  screen  discharge 
coefficient  much  less  than  0.1  which  when  inserted  into  Equation  18 
produced  high  values  of  the  screen  pressure  drop. 

The  technique  used  by  McDonnell  Douglas  to  calculate  pressure  drop 
across  surface  tension  screens  demonstrates  the  effect  of  viscosity  more 
directly.  The  current  correlation  used  by  McDonnell  Douglas  was  developed 
by  the  Win  Tec  Corporation  under  a NASA  contract  (Reference  7).  The 
correlation  is: 


APS  = apQ2  + bpQ  (21) 

2 

where,  AP$  = screen  pressure  drop,  Ib/in 

a,b  = constants  for  specific  screen  size  and  type 
p = liquid  density,  gm/ml 

p = liquid  viscosity,  centipoise 

Q = volumetric  flow  rate  per  unit  area,  9a^°ns  per  minute 

in^ 


Using  Equation  21  for  the  two  pressure  drop  calculation  examples 
previously  given  produces  the  results  shown  below: 


EXAMPLE  1:  44  micron  screen  at  1 g with  fuel  temperature  of  77°F 
Constants  a and  b are  obtained  from  Reference  8. 
a = 0.027 


b = 0.070 


n _ 10  lb-m  Y 1 

^ sec  A 9.0  lb-m 
gallon 


X 60 


sec 

min 


1 

160  in 


q _ 600  gallons  _ Q 41g7  gallons 
1400  min  sec^  min  in^ 

p = 1.07  gm/ml 

p = 24.8  centipoise 
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AP$  = 0.027  x 1.07  x (0.4167)2  + 0.070  x 24.8  x 0.4167 
AP$  = 0.00502  + 0.72339 
AP  = 0.72841  lb/in2 

5 

This  pressure  can  be  converted  to  inches  of  liquid  as  follows: 

. _ Ps  _ 0.72841  lb/in2 

s pg  0.03856  lb  x 1 

in3 

h$  = 18.89  inches 

EXAMPLE  2:  Fuel  Temperature  = 0°F 
Therefore,  p = 1.101  gm/ml 

u = 231  centi poise 

AP$  * 0.027  x 1.101  x (0.4167)2  + 0.070  x 231  x 0.4167 
AP$  = 0.00516  + 6.73804 
AP$  = 6.7432  lb/in2 

Converting  to  liquid  head, 

hs  ■ ffSnrr  ■ ,69-5  inches 

These  two  examples  demonstrate  the  effect  of  viscosity  on  screen  pressure 
drop.  Decreasing  the  temperature  from  77°F  to  0°F  increased  the  screen 
pressure  drop  by  nearly  an  order  of  magnitude.  Even  at  77°F  the  screen 
pressure  drop  obtained  by  this  technique  exceeds  the  total  permissible 
pressure  drop  (13.3  inches)  for  this  surface  tension  screen.  (See 
original  Example  1 calculation).  The  effect  of  the  other  parameters  in 
Equation  20  on  screen  pressure  drop  is  much  less  drastic  than  viscosity. 
The  constants  a and  b will  vary  with  specific  screen  type  but  will  vary 
only  by  a factor  of  about  two  for  the  constant  a and  about  four  for  the 
constant  b for  the  screens  used  in  this  study.  Viscosity  varies  by  a 
factor  of  about  nine  between  77°F  and  0°F  and  of  course  becomes  much 
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greater  as  the  fuel  temperature  is  decreased  toward  -65°F.  The  variance 
of  density  is  negligible  compared  to  the  variance  in  the  other  parameters. 
The  final  parameter  is  volumetric  flow  rate  per  unit  area.  This  parameter 
will  be  less  than  one  unless  a fairly  small  screen  is  used.  Therefore, 
this  parameter  has  much  less  effect  on  screen  pressure  drop  than  does 
viscosity  since  it  is  squared  in  the  first  term  of  Equation  21  and  is 
multiplied  by  viscosity  in  the  second  term  of  Equation  21.  Thus,  it  is 
seen  that  viscosity  is  the  primary  parameter  controlling  pressure  drop 
across  surface  tension  screens  when  one  is  considering  a high  viscosity 
fluid  such  as  RJ-5. 
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SECTION  III 
TRAP  TANK  SIZING 

1 . TRAP  TANK  GEOMETRY  AND  CALCULATION  PROCEDURES 

During  normal  positive  acceleration  flight  conditions  the  trap  tank 
containing  the  surface  tension  screens  is  completely  covered  and  filled 
with  fuel  so  that  fuel  movement  (sloshing)  in  the  main  tank  has  minimal 
effect  on  the  operation  of  the  surface  tension  units.  During  certain 
maneuvers  of  the  missile,  however,  negative  g conditions  will  exist  which 
will  uncover  the  trap  tank  for  short  periods  of  time.  The  essential 
requirement  for  the  trap  tank  is  to  supply  sufficient  fuel  for  the  length 
of  time  of  the  negative  g maneuver  so  that  the  engine  will  not  flame  out. 
An  analysis  of  the  retention  time  provided  by  trap  tanks  for  use  in  the 
Harpoon  missile  is  presented  in  Reference  3.  This  analytical  technique 
has  been  adapted  for  use  with  RJ-5  and  is  presented  in  the  following 
paragraphs. 

The  trap  tank  flow  model  used  for  the  retention  time  analysis  is 
shown  in  Figure  7.  The  total  volumetric  fuel  flow  rate  out  of  the  tank 
is  given  by 


^ = | + A 7206  (22) 

where  4r-  = total  volumetric  flow  rate,  ft^/sec 
m = engine  mass  flow  rate,  lb-m/sec 

3 

p = fuel  density,  lb-m/ft 

2 

A = total  effective  orifice  area,  ft 

2 

g = gravitational  acceleration,  32.17  ft/sec 

h = fluid  height  above  orifice,  ft 

G = number  of  negative  G's,  dimensionless 
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Assuming  that  the  fluid  volume  is  linear  with  fluid  height  gives, 

V = a + kh  (23) 

where  a and  k are  constants.  Differentiating  Equation  23  and  noting  that 
fluid  height  decreases  with  time  gives 


dv  _ . dh 
dt  " _K  dt 


Equating  Equations  22  and  24  gives 


(24) 


-k 


dh  _ m 
dt  ~ p 


+ A72g“hG 


(25) 


Integration  of  Equation  25  and  insertion  of  the  boundary  conditions  of 
h = hQ  at  t = 0 gives  the  following  expression  relating  fluid  height  with 
time:. 

(26) 

where  a = Av/^gG” 


t = ^ 

L 9 


1/? 

1 /o  1/9  / Qq  ^ ^ 

“(h°  _h  ’ * Q*  ,n(  ) 


The  initial  calculation  made  was  for  the  volume  of  fuel  left  in  the  trap 
tank  as  a function  of  fuel  height.  The  configuration  of  the  trap  tank 
considered  for  this  analysis  is  shown  in  Figures  7 and  8.  From  the 
schematic  shown  in  Figure  8,  the  following  geometrical  relationships  are 
derived: 


Cos 

a- 

n 

CD 

(27) 

Sin 

fl  - c/2  _ c 
s " R " 2R 

(28) 

S 

= 2R  e 

(29) 

H 

= R(l-Cos  e) 

(30) 

C 

= 2R  Sin  0 

(31) 
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Figure  7.  Trap  Tank  Flow  Model 
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Figure  8.  Trap  Tank  Schematic 
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where, 

R = Radius  of  main  tank,  inches 

H = height  of  liquid  in  trap  tank,  inches 

C = width  of  liquid  in  trap  tank,  inches 

6 = half-angle  formed  by  radii  drawn  from  main  tank  center  to 
the  intersection  of  liquid  level  with  trap  tank  wall,  radians 

The  volume  of  liquid  is  given  by  the  following  equation: 

VL I QU I D = Lr2(0  “ 1/2  Sin  20> 

of  ^LIQUID  = ? [RS-C(R-H)]  (32) 

L = Trap  tank  length,  inches 

Variation  of  fuel  volume  as  a function  of  fuel  height  was  obtained  for 
two  trap  tank  heights.  The  first  trap  tank  considered  was  a 9-inch  trap 
tank  height  with  the  main  fuel  tank  being  18  inches  in  diameter.  Trap 
tank  lengths  of  6,  9,  and  12  inches  were  used  to  obtain  three  trap  tank 
sizes.  Results  of  the  calculations  for  th.  first  trap  tank  height  are 
presented  in  Table  II  and  shown  graphical!  by  Figure  9.  The  data  points 
plotted  in  Figure  9 are  solutions  to  Equation  32.  These  data  are  then 
approximated  by  two  linear  segments  for  each  tank  size  so  that  slopes 
and  the  point  where  the  slopes  change  could  be  obtained  for  each  tank 
size.  The  linear  approximation  provides  the  constant  required  in 
Equation  23.  As  will  be  shown  later,  this  information  is  necessary  for 
obtaining  fuel  expulsion  times  from  the  trap  tanks. 

The  second  trap  tank  considered  had  a height  of  3.2  inches  with  a 
length  of  16  inches.  The  main  tank  diameter  was  maintained  at  18  inches. 
Table  III  presents  the  fuel  volume  as  a function  of  fuel  height  for  this 
tank.  Figure  10  shows  the  trap  tank  fuel  volume  as  a function  of  fuel 
height.  Again,  the  data  obtained  was  plotted  as  two  linear  segments  so 
that  fuel  expulsion  times  could  be  obtained. 
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TABLE  II  - FUEL  VOLUME  VS.  FUEL  HEIGHT  FOR  NINE-INCH  TRAP  TANKS 
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TABLE  III  - FUEL  VOLUME  VS.  FUEL  HEIGHT  FOR 


3.2  INCH  HIGH  TRAP  TANK 


FUEL  HEIGHT,  INCHES 

3.2 

3.0 

2.6 

2.2 

2.0 

1.3 

1.6 

1.4 
1.2 
1.0 
0.9 
0.8 
0.6 
0.4 


FUEL  VOLUME , INCHES 

489.5 
446.0 

352.6 

284.3 
2^7.3 
211.9 
178.2 

146.4 

116.6 

89.0 

76.1 
63.9 

41.6 

22.7 
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Figure  10.  Trap  Tank  Fuel  Volume  Versus  Fuel  Height  for  3.2  Inch  Tank 
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Equation  26  is  a general  expression  relating  fuel  height  in  the  trap 
tank  as  a function  of  time.  This  analysis  was  primarily  concerned  with 
the  total  time  required  to  empty  the  trap  tank  as  a function  of  engine 
mass  flow  rate.  An  expression  for  maximum  trap  tank  duration  can  be 
obtained  by  inserting  the  expression  for  linear  function  of  height  and 
volume  (Equation  23)  into  the  general  relationship  for  trap  tank  duration 
(Equation  26)  and  by  applying  the  boundary  conditions  that  t=t  when 
h = 0.  The  resulting  equation  obtained  is: 


max 


= - [k, 

a L 1 


(h0’«  - h,>'2)  ♦ k2  h,1'2] 


f 


k^ln 


(Qe  + h^a) 
♦ ho'/2“) 


+ k2ln 


( Qe  > 

(Qe  + h1/?a) 


(33) 


where,  k^  and  k^  ° slopes  of  lines  from  either  Figure  9 or  10 
hQ  = initial  fluid  height  in  tank,  inches 
h^  = fluid  height  at  point  of  slope  change,  inches 


Equation  33  was  programmed  for  solution  by  the  Hewlett-Packard  Model 

9100  Calculator.  For  each  run,  all  the  terms  on  the  right  hand  side  of 

Equation  33  are  input  as  constants  including  an  initial  value  of  engine 

mass  flow  rate  (m).  The  program  then  increments  m by  0.5  and  calculates 

t for  each  mass  flow  rate  until  the  machine  is  stopped  by  the  operator, 
max 

Complete  details  of  the  computer  program  are  described  in  the  Appendix. 


2.  RESULTS  FOR  9 INCH  AND  3.2  INCH  TRAP  TANKS 

Typical  results  are  presented  in  Figures  11  through  17.  Figure  11 

shows  the  effect  of  orifice  vent  area  for  the  9-inch  high  by  9-inch  long 

trap  tank  for  a sustained  load  of  negative  ten  g's.  Orifice  vent  areas 
2 2 2 

of  0.15  in  , 0.20  in  , and  0.25  in  were  considered  over  a fuel  flow 
rate  range  from  1.0  to  16  lb/sec.  The  retention  times  for  these  three 
orifice  vent  areas  were  significantly  different  from  one  another  only 
at  flow  rates  below  3.0  Ib/sec.  Above  this  flow  rate  the  difference  in 
retention  time  was  less  than  one  second.  The  other  factor  to  be 
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Figure  11.  Tank  Emptying  Time,  9 Inch  Height,  9 Inch  Length, 
10  6 1 s , Various  Vent  Areas 
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considered  is  the  total  retention  time  obtained,  regardless  of  orifice 
vent  area  considered,  as  a function  of  fuel  flow  rate  to  the  engine. 
Figure  11  shows  that  the  retention  times  for  all  three  orifice  vent  areas 
decreases  rapidly  with  increasing  fuel  flow  rate.  At  1.0  Ib/sec  engine 
flow  rate  the  retention  times  range  from  16.2  to  21.5  seconds.  At 
10  lb/sec  engine  flow  rate  the  retention  time  was  less  than  4 seconds  for 
all  orifice  vent  areas. 

The  effect  of  acceleration  (g)  level  on  retention  capability  is  shown 

in  Figure  12.  Acceleration  levels  of  1,  5,  and  15  g's  were  evaluated 

with  the  nine  inch  high  by  nine  inch  long  trap  tank  using  an  orifice 
2 

vent  area  of  0.20  in  . Again,  as  with  the  orifice  vent  area  parameter, 
g level  has  a significant  effect  only  at  low  flow  rates.  At  flow  rates 
above  4 lb/sec  the  difference  in  retention  times  for  the  1 g and  the 
15  g cases  is  less  than  2 seconds.  For  the  1 g case,  retention  time 
decreases  rapidly  from  18  seconds  at  2 Ib/sec  to  less  than  4 seconds 
above  11  Ib/sec  flow  rate.  For  the  15  g case,  retention  time  at  2 lb/sec 
is  only  11.5  seconds  and  rapidly  decreases  to  under  4 seconds  above  fuel 
flow  rates  of  9 Ib/sec. 

The  effect  of  trap  tank  size  was  examined  by  first  decreasing  tank 

height  to  6 inches  and  maintaining  the  9-inch  tank  length.  Figure  13 

shows  the  result  of  this  change  on  retention  time  for  the  three  orifice 

vent  areas  previously  discussed  and  also  at  a 10  g acceleration  level. 

For  the  0.15  in  area  at  1 Ib/sec  flow  rate  the  retention  time  is 

14  seconds  versus  21.5  seconds  for  the  9 inch  high  trap  tank.  At 

2 

10  Ib/sec  the  retention  time  for  the  0.15  in  case  is  about  2.4  seconds. 

2 2 

Similar  decreases  in  retention  time  for  the  0.20  in  and  0.25  in  orifice 
vent  areas  were  obtained  for  the  6 inch  high  trap  tank  versus  the  9 inch 
high  trap  tank. 
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Figure  12.  Tank  Emptying  Time,  9 Inch  Height,  9 Inch  Length, 

o 

0.20  In  Vent  Area,  Various  G Levels 
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Figure  13.  Tank  Emptying  Time,  6 Inch  Height,  9 Inch  Length, 
10  G's,  Various  Vent  Areas 
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Figure  14  presents  the  effect  of  g level  on  trap  tank  duration  for  the 
6 inch  high  by  9 inch  long  trap  tank.  For  this  tank  the  effect  of  g level 
becomes  important  at  fuel  flow  rates  below  3 Ib/sec.  At  this  flow  rate 
the  difference  in  retention  times  between  the  1 g case  and  the  15  g case 
is  less  than  two  seconds.  Retention  times  at  1 lb/sec  fuel  flow  rate 
ranged  from  19  seconds  for  the  1 g case  to  11  seconds  for  the  15  g case. 

At  10  lb/sec,  retention  time  for  both  g levels  was  about  2.5  seconds. 

The  second  variation  made  to  change  trap  tank  size  was  to  increase 
length  while  maintaining  the  9 inch  height  dimension.  The  same  type  of 
curves  as  before  were  obtained  for  the  three  orifice  vent  areas  and  for 

2 

the  three  acceleration  levels  at  a constant  orifice  vent  area  of  0.20  in  . 
The  primary  effect  was  tc  increase  retention  time  over  the  9 inch  length 
tank  by  a factor  equal  to  the  change  in  tank  volumes.  Data  for  the  12  inch 
length,  9 inch  high  trap  tank  are  plotted  in  Figures  15  and  16. 

Figure  17  presents  retention  time  data  for  all  three  tanks  having  the 

common  9 inch  height  dimension.  Tank  lengths  considered  were  6,  9 and 

12  inches.  These  calculations  were  made  for  the  10  g level  using  an 

2 

orifice  vent  area  of  0.20  in  . The  retention  time  versus  flow  rate  curves 
obtained  were  the  same  type  as  before  and  the  difference  in  retention 
times  was  a direct  function  of  the  tank  length.  Thus,  the  retention  time 
for  the  12  inch  length  tank  at  any  fuel  flow  rate  is  very  nearly  twice  the 
retention  time  for  the  6 inch  length  tank.  Rete*  iion  time  for  the  12  inch 
length  trap  tank  decreased  from  32  seconds  at  a fuel  flow  rate  of  0.5  lb/sec 
to  about  5 seconds  at  10  lb/see  flow  rate  whereas  for  the  6 inch  tank 
length  the  retention  time  decreased  from  16.6  seconds  at  0.5  lb/sec  flow 
rate  down  to  about  2.6  seconds  at  10  lb/sec  flow  rate. 
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Figure  14. 


Tank  Emptying  Time, 
0.20  In^  Vent  Area, 


6 Inch  Height,  9 Inch  Length, 
Various  6 Levels 
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Figure  ’5.  Tank  Emptying  Time,  9 Inch  Height,  12  Inch  Lengi.h, 
10  6' s Various  Vent  Areas 
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Figure  16.  Tank  Emptying  Tin^ , 9 'inch  Height,  12  Inch  Length, 
0.20  In2  Vent  Area,  Various  G Levels 


\ 
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The  dimensions  and  volumes  of  the  four  trap  tanks  considered  so  far 
by  this  analysis  are  as  follows: 


TABLE  IV  — TRAP  TANK  VOLUMES 
Trap  Tank  dimensions  Volume , In^ 


9-inch  height,  6 -inch 
length 

763 

9-inch  height,  9-inc'n 
length 

11U5 

9rinch  height,  12-inch 
length 

1527 

6-inch  height,  9-inch 
length 

670 

As  shown  by  Figures  11  through  17,  the  retention  times  available 

with  the  four  tanks  listed  in  Table  IV  vary  from  over  40  seconds  at  fuel 

flow  rates  of  1.0  lb/sec  and  less  to  under  three  seconds  for  flow  rates 

of  10  lbs/sec  and  over.  Specific  retention  time  available  from  a given 

trap  tank  will  also  depend  on  G level  and  orifice  vent  area  as  previously 

discussed.  Fuel  volume  versus  fuel  height  for  the  3.2  inch  high  trap 

tank  are  presented  in  Table  III  and  shown  graphically  in  Figure  10. 

3 

Maximum  fuel  volume  for  this  trap  tank  is  490  in  . Maximum  trap  tank 

retention  times  were  calculated  for  this  smaller  tank  using  Equation  33. 

Figures  18  and  19  show  the  effect  of  orifice  vent  area  at  the  10  g 

2 2 

acceleration  level.  Orifice  vent  areas  used  were  0.10  in  , 0.15  in  , 

2 2 

0.20  in  , and  0.30  in  . Again,  the  effect  of  orifice  vent  area  was 
significant  only  at  flow  rates  below  2 lb/sec.  At  2 Ib/sec  the  difference 
in  retention  time  for  the  0.10  in  case  and  0.30  in  case  was  less  than 
two  seconds.  As  the  flow  rate  was  decreased  toward  zero  the  retention 
time  difference  for  the  two  vent  areas  became  significant  whereas 
increasing  flow  rate  above  2 lb/sec  decreased  the  difference  to  an 
insignificant  level. 
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Figure  18.  Tink  Emptying  Time,  3.2  Inch  Height,  16  Inch  Length, 
10  G's,  0.10  and  0.20  In2  Vent  Areas 
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AFAPL-TR-74-29 


Figure  20  shows  the  effect  of  G level  on  retention  time  for  the 
3.2  inch  high  trap  tank.  Orifice  vent  area  used  was  0.20  in  . Acceler- 
ation levels  of  1,  5,  and  20  g's  were  used  to  calculate  retention  time 
as  a function  of  fuel  flow  rate  to  the  engine.  Once  again,  it  is  seen 
that  acceleration  level  is  significant  only  at  relatively  low  mass  flow 
rates.  At  4 Ib/sec  the  difference  in  retention  times  between  the  1 g and 
the  20  g cases  is  less  than  one  second  but  does  increase  to  an  appreciable 
difference  at  flow  rates  below  1.5  lb/sec.  The  retention  time  available 
with  the  3.2  inch  high  trap  tank  vary  from  about  14  seconds  with  the 
0.20  inch  orifice  vent  area  configuration  at  1.0  lb/sec  flow  rate  and 
1 g to  less  than  two  seconds  for  all  the  conditions  evaluated  at  flow 
rates  above  9.0  lb/sec.  These  retention  times  are  than  about  one-third 
of  the  9 inch  high  trap  tank  values  at  1.0  lb/sec  flow  rate  and  nearly 
the  same  values  for  flow  rates  above  10  lb/sec.  These  results  can  then 
be  used  to  make  a first  estimate  for  trap  tank  size  for  a specific 
application.  The  critical  factor  is  the  retention  time  required  at  low 
flow  rates  since  there  exists  a three-fold  factor  in  retention  times  at 

1.0  lb/sec  for  the  two  basic  trap  tank  configurations  considered.  If 
retention  times  of  less  than  14  seconds  at  low  flow  rates  are  adequate, 
the  3.2  inch  high  trap  tank  would  be  the  logical  choice  since  this  places 
less  restrictions  on  the  surface  tension  screen  sizes,  fuel  temperatures, 
and  acceleration  levels  which  can  be  accommodated.  At  flow  rates  above 

5.0  Ib/sec  the  trap  tank  size  becomes  less  critical  with  regard  to 
retention  time  since  the  absolute  values  of  the  retention  times  is  less 
than  ten  seconds  and  the  maximum  difference  in  retention  times  for  the 
two  primary  configurations  is  less  than  6 seconds. 

The  trap  tank  designs  considered  so  far  by  this  analysis  have  assumed 
that  in  the  operational  situation  the  axial  g's  are  always  unidirectional 
so  that  the  trap  tank  can  be  located  at  one  end  of  the  main  fuel  tank. 

This  would  normally  be  the  case  for  the  climb  and  cruise  segments  of  a 
missile  trajectory  where  the  acceleration  force  would  keep  the  fuel  at 
either  the  aft  end  and/or  bottom  of  the  main  fuel  tank.  If  missile 
maneuvering  is  required  near  the  end  of  a mission  when  the  main  fuel 
tank  is  nearly  empty  the  maneuvers  would  need  to  be  restricted  to  those 
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with  a positive  acceleration  force  in  the  direction  of  the  trap  tank. 

This  restriction  comes  about  due  to  the  hydrostatic  head  limitation  and 
configuration  of  the  trap  tank.  Assuming  that  the  maximum  hydrostatic 
head  which  can  be  sustained  is  3.2  inches  of  RJ-5,  the  location  of  the 
surface  tension  screen(s)  inside  the  trap  tank  must  be  configured  so  that 
the  3.2  inch  restriction  is  not  exceeded  during  any  missile  maneuver. 

If  this  limit  is  exceeded,  pressurant  gas  ingestion  will  occur  the 
moment  any  portion  cf  the  surface  tension  screen  becomes  uncovered  from 
the  fuel . 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

Surface  tension  devices  have  been  employed  in  a number  of  missiles 
and  space  vehicles  to  orient  the  fuel  and  to  insure  fuel  flow  during 
negative  g maneuvers  (Reference  I).  Most  of  the  applications  to  date 
have  been  for  small  (less  than  1)  negative  g maneuvers  and  the  fuel 
(or  propellant)  has  been  a low  viscosity  liquid.  The  analysis  performed 
in  this  study  considered  much  higher  negative  g levels  than  did  previous 
applications  of  surface  tension  screens  and  also  considered  a fuel  which 
is  several  orders  of  magnitude  more  viscous  (at  low  temperatures)  than 
the  JP-type  fuels  or  rocket  propellants  used  in  previous  fuel  systems. 

This  extreme  combination  of  a viscous  fluid  at  a high  negative  g level 
make  the  use  of  surface  tension  devices  impractical  if  not  impossible. 
However,  at  low  fuel  flow  rates  and  moderate  (up  to  3)  negative  g levels, 
surface  tension  devices  may  be  applicable  even  with  the  viscous  RJ-5 
fuel.  This  would  depend  to  a great  extent  on  keeping  the  fuel  temperature 
at  a moderately  high  level,  perhaps  above  40°F  as  an  example.  This 
temperature  restriction  can  be  relieved  by  using  a less  viscous  fuel  such 
as  the  JP-9  type  fuel  currently  being  developed  or  by  heating  the  RJ-5 
fuel.  If  further  analysis  of  the  MPM  reveals  that  operation  in  a low 
negative  g environment  is  feasible,  the  utilization  of  surface  tension 
screens  may  prove  to  be  practicable. 
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APPENDIX 

CALCULATION  TECHNIQUE  FOR  TRAP  TANK  SIZING 

The  basic  assumption  used  for  the  trap  tank  sizing  analysis  was  that 
the  fluid  height  remaining  in  the  tank  as  a function  of  time  is 
approximated  by  a linear  relationship  of  the  form: 


V = a + kh 


(34) 


Differentiation  of  this  equation  and  combination  of  the  resulting 
differential  equation  with  the  expression  for  total  volumetric  flow  rate 
out  of  the  tank. 


dV 

dt 


= + AV2ghG 


(35) 


produces  the  following  expression  for  expulsion  time: 


-hl1/2)  + 2Me 
a 


In 


5.  ♦ “hoV7 


(36) 


To  obtain  more  nearly  exact  volumes  as  a function  of  fluid  height, 
Equation  34  is  evaluated  over  two  segments  as  indicated  by  Figure  21. 
Total  expulsion  time  then  is, 


total  = tl 


+ t„ 


(37) 
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FLUID  HEIGHT,  INCHES 

Figure  21.  Fluid  Volume  Vs  Fluid  Height 


The  total  expulsion  time  for  each  tank  geometry  is  obtained  by 
determining  the  two  slopes  (k^  and  k2)  of  the  linear  approximation  lines 
shown  by  Figure  21  while  keeping  the  main  tank  diameter  constant  and 
using  Equation  36  as  follows: 


*,  • £ - *,,/2) 


2kl  Qe  In 


a 


zh.  (h1,/2-n2,/2)*^i,n 


(38) 


(39) 


Since  h2  approaches  zero  as  the  fluid  volume  approaches  zero,  Equation 
39  reduces  to: 


t2  = 


2k9 

— - h, 
a 1 


1/2  + 2k2Qe 


In 


Qe  + ah1 


1/2 


(40) 
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From  Equation  37, 


Hotal  ^max  *"1  + *2 


(41) 


Substituting  Equations  38  and  40  into  Equation  41  and  simplifying  gives: 


fmox  = •§•  [K,(h0*  “hi  *)  +K2hi1E]  +^[K|  ln(~aT"L!i)  +K2 ln  ( — I')] 

aL  J CT  L 'Q.+CrtvJ'  '0.+ah.I/J 


Qa+ah|^> 


O.+ahn 


'Q.+ah.? 


(42) 


where  a = A -\/2gc6  = A ./772.8G 


All  terms  on  right  hand  side  of  Equation  42  can  be  input  as  constants 
for  one  run  including  initial  value  of  m.  m is  incremented  for  each  time 
calculation  by  adding  0.5.  Program  will  continue  to  calculate  t „ for 
each  m until  machine  is  stopped  by  operator. 

Program  1 and  Program  2 are  recorded  on  a magnetic  card  for  reuse. 

Before  executing  program,  store  values  in  registers: 

A + 201 

G + 202 

hn  + 203 
°1 

k1  - 204 
k2  -»•  205 
h1  206 
p -v  207 

m + 200  (initial  value) 
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PROCEDURE 

1.  Store  constants  in  registers  as  on  previous  page 

2.  Put  Program  1 and  Program  2 in  extended  memory 

3.  Turn  on  printer  with  "X"  key  depressed 

4.  Press  1,  FMT,  GO  TO,  CONTINUE 

5.  Program  will  continue  until  operator  presses  STOP 

ADDITIONAL  CASES 

1.  Enter  any  new  constants  in  registers 
NOTE:  Reset  initial  value  of  m 


2.  Start  with  instruction  4 above 

OUTPUT  FORMAT 
SEE  RIGHT 


A = 

0.20 

G = 

1.0 

1)  The  program  first 

V 

9.0 

outputs  the  initial 

constants  stored  in  the 

hl  = 

3.0 

registers. 

kl  = 

200 

k2 

100 

P = 

0.46 

m = 

0.5 

2)  Then  values  of  m 

t = 

130.323- 

and  its  corresponding 

1.0 

t are  printed  as  the 

m = 

program  continues  to 

t = 

116.557 

increment  m and  solve 

. 

1.5 

for  t. 

m = 

106.121 

t = 

1 

1 
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-nut  Program  1 (cont.) 
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Title  ProSram  2 (Subroutine) 
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